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Models in which dark matter particles can scatter into a slightly heavier state which promptly 
decays to the lighter state and a photon (known as eXciting Dark Matter, or XDM) have been shown 
to be capable of generating the 3.55 keV line observed from galaxy clusters, while suppressing the flux 
of such a line from smaller halos, including dwarf galaxies. In most of the XDM models discussed in 
the literature, this up-scattering is mediated by a new light particle, and dark matter annihilations 
proceed into pairs of this same light state. In these models, the dark matter and mediator effectively 
reside within a hidden sector, without sizable couplings to the Standard Model. In this paper, we 
explore a model of XDM that does not include a hidden sector. Instead, the dark matter both 
up-scatters and annihilates through the near resonant exchange of a 0(10^) GeV pseudoscalar with 
large Yukawa couplings to the dark matter and smaller, but non-neglibile, couplings to Standard 
Model fermions. The dark matter and the mediator are each mixtures of Standard Model singlets 
and SU{2)w doublets. We identify parameter space in which this model can simultaneously generate 
the 3.55 keV line and the gamma-ray excess observed from the Galactic Center, without conflicting 
with constraints from colliders, direct detection experiments, or observations of dwarf galaxies. 

PACS numbers: 95.35.+d, 95.85.Pw; FERMILAB-PUB-15-009-A 


I. INTRODUCTION 

The nature of dark matter remains one of the most 
elusive and longstanding problems in physics today. As 
a consequence, much attention has been given to obser¬ 
vational anomalies that can be plausibly interpreted in 
terms of dark matter interactions. One such signal is 
an approximately 3.55 keV X-ray line that has been ob¬ 
served from a number of galaxy clusters, as well as from 
the nearby Andromeda Galaxy. 

The first reported evidence for the 3.55 keV line was 
found in data from the XMM-Newton satellite, from the 
directions of a stacked sample of 73 low redshift galaxy 
clusters [1]. Shortly thereafter, a similar line was re¬ 
ported from the directions of the Perseus Cluster and the 
Andromeda Galaxy [2]. A study of XMM-Newton data 
also suggests the existence of a 3.55 keV line from the di¬ 
rection of the Milky Way’s center [3] (see also, however. 
Ref [4]). More recently, the line was identified within 
Suzaku data from the Perseus Cluster [5]. 

A number of interpretations for these observations 
have been proposed. On the one hand, it has been sug¬ 
gested that atomic transitions (such as those associated 
with the chlorine or potassium ions, Cl-XVII and K- 
XVIII, for example [6]) might be responsible for the line, 
although the viability of this explanation is currently un¬ 
clear [7-9]. Alternatively, decaying dark matter particles 
could generate such an X-ray line. Particularly well mo¬ 
tivated is dark matter in the form of an approximately 
7 keV sterile neutrino, which decays through a loop to 
a photon and an active neutrino. If one assumes that 
all of the dark matter consists of 7 keV sterile neutri¬ 
nos, the observed X-ray line flux implies a mixing angle 


of sin^(20) ~ 7 X 10“^^. With such a small degree of 
mixing, however, the standard Dodelson-Widrow mecha¬ 
nism of production via the collision-dominated oscillation 
conversion of thermal active neutrinos [10] leads to an 
abundance of sterile neutrinos that corresponds to only a 
few percent of the total dark matter density, thus requir¬ 
ing additional resonant or otherwise enhanced production 
mechanisms. Alternatively, sterile neutrinos with a larger 
mixing angle of sin^(20) ~ 3 x 10“^° could naturally con¬ 
stitute roughly 10% of the dark matter abundance, and 
decay at a rate that is sufficient to generate the observed 
line flux. 

Interpretations of the X-ray line in terms of decaying 
dark matter are in considerable tension, however, with 
studies of galaxies using Ghandra and XMM-Newton 
data [11] and dwarf spheroidal galaxies using XMM- 
Newton data [12], which do not detect a line at the 
level predicted by decaying dark matter scenarios. One 
way to potentially reconcile the intensity of the line ob¬ 
served from clusters with the null results from dwarfs and 
other smaller systems is to consider the class of scenarios 
known as eXciting Dark Matter (XDM) [13-16]. In such 
models, the collisions of dark matter particles can cause 
them to up-scatter into an excited state, XiXi -7- X 2 X 2 or 
XiXi ~7- XiX 2 - For a mass splitting of ~ '^xi — 
keV, the subsequent decays of the slightly heavier state 
can generate a 3.55 keV photon, X 2 —> Xi7- Critical 
to the problem at hand are the kinematics of the XDM 
scenario, which introduce a velocity threshold for up- 
scattering, suppressing the X-ray flux from dwarf galax¬ 
ies (and, to a lesser extent, from larger galaxies) [15, 17]. 
Within the paradigm of XDM, the observations of clus¬ 
ters, galaxies, and dwarf galaxies can be mutually con- 
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sistent for dark matter masses between approximately 40 
GeV and 10 TeV [17], covering the mass range generally 
associated with conventional WIMPs. 

If up-scattering WIMPs are responsible for the 3.55 
keV line, one might also imagine that the same dark 
matter species could generate the excess of GeV-scale 
gamma-rays observed from the region surrounding the 
Galactic Center [18-27]. This signal, identified within 
data from the Fermi Gamma-Ray Space Telescope, ex¬ 
hibits a spectrum and morphology that are in good agree¬ 
ment with that anticipated from dark matter annihila¬ 
tions. This data has been explored by several groups 
independently, including recently the Fermi Collabora¬ 
tion [28]. Assuming annihilations to bb, for example, dark 
matter particles with a mass of ~ 35-65 GeV and a 
cross section of (av) ~ 10“^® cm®/s provide a good fit to 
the observed excess [29]. 

The primary challenge in developing a viable XDM 
model for the 3.55 keV line is that the up-scattering rate 
must be very high, several orders of magnitude larger 
than the annihilation rate. One way to realize this is to 
consider dark matter that scatters through a light media¬ 
tor and annihilates into pairs of the same mediator. This 
naturally leads to an up-scattering rate that is enhanced 
by a factor of ~ {my,/MMed)‘^ relative to the annihila¬ 
tion rate. As this phenomenology can be realized with¬ 
out the dark matter or mediator possessing any sizable 
couplings to the Standard Model (SM), these scenarios 
are sometimes called “hidden sector” models. Examples 
of such proposals include models with a massive vector 
(hidden photon) or a massive scalar (hidden Higgs) that 
couples directly to the dark matter, but interacts with 
the SM only through a very small degree of kinetic or 
mass mixing. As a result, the dark sector and SM are 
effectively sequestered from one another. As this class 
of possibilities has been explored previously in some de¬ 
tail [15, 16, 30-32], we do not consider it here. Instead, 
we explore models in which the dark matter annihilates 
directly into SM fermions (for an earlier investigation in 
this direction, see Ref. [33]). By introducing a resonant 
mediator with a hierarchy of couplings ((/dark >> 5 sm)) it 
is possible to accomplish similar phenomenology without 
a light mediator. We identify such a model that can si¬ 
multaneously explain the 3.55 keV line from Galaxy Glus- 
ters and the Galactic Center gamma-ray excess. We find 
viable parameter space in our model that is consistent 
with all current collider, direct detection, and indirect 
detection constraints. 


II. THE KINEMATICS OF EXCITING DARK 
MATTER 


If the 3.55 keV signal is due to dark matter, the model 
responsible needs to address why this signal is not seen 
from dwarf galaxies (and, to a lesser extent, from larger 
galaxies). As the up-scattering rate in the XDM scenario 
depends strongly on the dark matter velocity dispersion 


in such systems, this framework provides a simple mech¬ 
anism to suppress the line flux predicted from smaller 
halos. 

The velocity averaged cross section for up-scattering is 
given by: 


{av) = aova, (1) 

where the normalization, (jg, is taken to be a free param¬ 
eter and 7 accounts for the effect of the threshold velocity 
on the up-scattering rate: 


1 = - I Q{v-vt)^. ( 2 ) 

The quantity Vt is the threshold velocity, given by: 


Vt 



( 3 ) 


where (taken to be ~ 3.55 keV) is the mass splitting 
between X 2 and xi and A = 1 (2) for up-scattering to 
XlX2 (X2X2)- 

In the limit of 5m^ = 0 (and Vt = 0), the standard 
{av) = (Tgu is recovered. For larger mass splittings, how¬ 
ever, the up-scattering rate and corresponding line flux 
will be suppressed in smaller systems, where typical ve¬ 
locities are lower. 

To obtain up-scattering rates in dwarfs, galaxies, and 
clusters that are each compatible with the reported ob¬ 
servations, Ref. [17] finds that a threshold velocity of 
Vt — 2Q — 245 km/s is required (at the < 3 level, 
and assuming that the excited state decays promptly). 
Combining this with Eq. 3 (where N = 2), this implies 
~ 40 GeV-10 TeV. Annihilating dark matter parti¬ 
cles near the low end of this mass range are also well 
suited to account for the Galactic Center gamma-ray ex¬ 
cess. 


III. MODEL BUILDING 
A. Up-scattering 

There are two classes of scenarios in which an excited 
state could be presently decaying in order to generate the 
observed 3.55 keV line. First, if the excited state has a 
lifetime on the order of the age of the Universe or longer, 
a population of such particles could have been produced 
in the early universe. Primordial excitations, however, 
do not lead to a relative suppression in dwarf galaxies, 
and thus suffer from the same challenges in explaining the 
3.55 keV line as ordinary decaying dark matter. Alterna¬ 
tively, if the excited state is short lived (millions of years 
or less) collisions between dark matter particles must lead 
to an up-scattering rate that is sufficient to perpetually 
populate these excitations in galaxy clusters. It is this 
second case that we consider here. 
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FIG. 1. The dominant Feynman diagram for up-scattering in 
our model. 


In order for XDM to generate the flux of 3.55 keV 
photons observed from galaxy clusters, very large cross 
sections for up-scattering are required, in the approx¬ 
imate range of <Jv{xiXi X 2 X 2 ) ~ GeV) x 

10“^® cm^/s. In addition to being very large in and 
of itself, this value for the up-scattering cross section 
is several orders of magnitude larger than the annihila¬ 
tion cross section needed to generate the Galactic Genter 
gamma-ray excess, or to obtain a thermal relic abundance 
in agreement with the measured dark matter density. 

In light of this, it is interesting to consider the up¬ 
per limit imposed on dark matter scattering from the 
point of view of perturbativity and unitarity. In this 
paper, we will focus on up-scattering through a resonant 
s-channel pseudoscalar, a (see Fig. 1). We will further as¬ 
sume that the dark matter and its excited state, Xi, 2 j are 
each Majorana fermions with nearly degenerate masses, 
(collectively constituting a pseudo- 
Dirac fermion). The scalar (</ = 0) bilinear involved 
in this interaction, x^l^Xi being even under charge, 
C = (—1)^''’'^, and odd under parity, V = (—1)^"'’^, im¬ 
plies that this operator only acts on incoming dark mat¬ 
ter pairs with zero spin and orbital angular momentum, 
J = S = L = 0. As a result, scattering through an s- 
channel pseudoscalar is purely s-wave, and the unitarity 
bound (see e.g. Ref. [34]) on up-scattering is given by 
(assuming X 2 is self-conjugate): 


27r 

av < —^ 
m^v 


/50 GeV 

V 


0.003 


V 


xlO ^^cm^/s, (4) 


where v ~ 0.003 is the typical dark matter relative ve¬ 
locity in a galaxy cluster. For comparison, note that 
this upper limit is much stronger than that derived 
from self-scattering in objects such as the Bullet Clus¬ 
ter, <JV < ( 50 ^) (oib) 10-14 crnVs [35]. 

A more explicit bound on the couplings of the theory 
arises if one parametrizes the Lagrangian responsible for 
up-scattering as follows: 

£ D Xhaxiix^Xi + A22aX2*7®X2 + A?2aXD7®X2- (5) 


If we assume that xi and X 2 couple to the pseudoscalar, 
a, with approximately equal strength, we can use Eq. 4 


to deduce: 


Ai 1 22 ^ y/27r 


{s-mjf + mlTl 
s (s — 4m^) 


1/4 


( 6 ) 


where nia and Fq are the mass and width of a, respec¬ 
tively. Therefore, perturbative unitarity of the theory 
in the non-relativistic and ultra-relativistic regimes re¬ 
quires: 


^11,22 


< 


20 X (d/O.l)^/^ (0.003/v)^/% 
2.5, 


y/s ~ 2m^ 
~ 00, 


( 7 ) 

where d = jl — (ma/2m^) ] <C 1. In the first line of 
Eq. 7, we have assumed that the width Fa is sufficiently 
small such that Fa/ma ^ S. We will show later that 
these conditions will be satisfied within the most viable 
parameter space for generating a large cross section for 
up-scattering. We see that in the non-relativistic regime, 
the upper limit on 22 from perturbative unitarity is 
weaker than one generically expects from perturbativity 
of the theory, AJ^ 22 ^ ^tt. Throughout our analysis, we 
will consider Yukawa couplings as large as 22 ~ ^ in 
the non-relativistic regime. If xi ,2 only couple to a, then 
large values of A“i 22 ''vill contribute positively to its beta 
function and cause A“]^ 22 to grow rapidly at higher en¬ 
ergies. By considering such large values of this coupling, 
we implicitly require that new physics (such as couplings 
to new gauge bosons) come in at higher energies in order 
to stabilize A“]^ 22 ~ 2.5 in the high-energy limit. 

In the low velocity limit, the cross section for XiXi to 
become excited to X 2 X 2 is given by: 


crw(xiXi X 2 X 2 ) « 


X 


(X-7)''" 

2m^ [A“iA22]^ 
7r[(4TO2 - m2)2 -h mlTl] ’ 


( 8 ) 


where Vt is as defined in Eq. 3. The up-scattering of XiXi 
into X 1 X 2 is subdominant, but for a somewhat subtle rea¬ 
son. As we will see later in this paper, after diagonalizing 
into mass eigenstates, one of the fields xi or X 2 gener¬ 
ally has a mass term with the “wrong sign”, requiring 
a transformation x ij^X- This leaves the axi*7^Xi 
and ax 2 il^X 2 interaction terms as written in Eq. 5, but 
changes the mixed term into 0 x 1 X 2 , resulting in the sup¬ 
pression of the corresponding up-scattering cross section 
by an additional factor of — Vt)- 


B. Annihilation and Coannihilation 

In the previous subsection, we showed that the very 
large up-scattering rates required for the 3.55 keV line 
can be generated through the resonant exchange of a 
pseudoscalar, a, but at the cost of introducing 0 ( 1 ) 
Yukawa couplings into the dark matter sector. If xi is to 
be populated thermally in the early universe, however. 
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FIG. 2 . The dominant Feynman diagram for the annihilation 
or coannihilation of XiXi, X1X2, or X2X2- 


it must also have non-zero couplings to the SM. Further¬ 
more, if the couplings of a to the SM are comparable to its 
couplings to the dark matter, then the annihilation cross 
section during freeze-out will be many orders of magni¬ 
tude larger than that needed to generate a thermal relic 
abundance of xi consistent with ~ 0.12. Instead, 
there must be a large hierarchy between the couplings of 
a with the dark matter and with the SM. 

One way to generate a very small coupling between 
the a and SM fermions is through mass-mixing with 
the heavy pseudoscalar in a two-Higgs doublet model 
(2HDM). Such a scenario has been discussed previously 
within the context of the Galactic Center gamma-ray ex¬ 
cess [36]. The idea is to introduce a scalar potential in¬ 
volving a parity-odd singlet pseudoscalar, uq, along with 
a second Higgs doublet in the framework of a Type-II 
2HDM. The two Higgs doublets, each with hypercharge 
of -1-1/2, are denoted as and the corresponding 

pseudoscalar of the 2HDM sector is written as Aq. After 
Gq and Aq mix, the light and heavy mass eigenstates of 
the CP-odd sector will be written as a and A, respec¬ 
tively. 

The Higgs portal between the dark matter and the 
SM emerges from the trilinear interaction of the scalar 
potential involving oq, and More specihcally, the 
terms of the scalar potential relevant for the annihilation 
and coannihiation of xi ,2 are given by: 


vacuum expectation values. Once electroweak symmetry 
breaking is induced by (Hd^u) = Vd,ulV^, Hd and Hu can 
be written in terms of the scalar mass eigenstates of the 
2HDM potential: 


= -V2 spH+ + V2 C0G+ \ 

^/2\Vd-Sah + CaH-iS0Ao+iC/3GJ ’ 

h. = A( + i , (10) 

\>Jd + Cah SaH ICpAf) -\- ISfdGJ 

where h, H are the light and heavy CP-even Higgs 
bosons, the charged Higgs, G and G^ the neutral and 
charged Goldstones bosons, and Aq is the pseudoscalar 
of the 2HDM sector, and sp are the cosine and sine 
of /?, defined by tan /3 = Vu/vd and \/v^ + Vu = v = 246 
GeV, and Ca and Sa are the cosine and sine of the mass 
mixing angle of the CP-even scalars, a. We will remove 
all dependence on a by choosing to work in the align¬ 
ment limit throughout, where sin(/3 — a) = 1 and the h 
couplings are SM-like. We will further assume that the 
masses of Aq, H, and are decoupled, with values at 
a scale around 1 TeV. 

Mass mixing between the pseudoscalars ag and Aq is 
induced by the coupling Ba in Eq. 9. The mass-squared 
matrix of the CP-odd sector in the (oq, Aq) basis is writ¬ 
ten as: 


-^CP-odd 


-BaV 

-BaV 


( 11 ) 


where toaq is the mass of the the pseudoscalar Aq in the 
2HDM potential. Diagonalizing leads to the 

mass eigenstates a and A such that: 


ao 

Aq 


cos 0 sin 9 


2 ^ 

A = 7 : 

a,A 2 


— sin 9 cos 9J A 


( 12 ) 




COS 0 = —= I 1 




4H>2 

1/2 


{mX - ml^ 


4H?r;2 


Hscaiar A V 2 HDM + + (iBaGgHlHu + h.c.^ , (9) 

where V 2 HDM is the most general CP-conserving 2HDM 
potential corresponding to a Type-H 2HDM, and Ba is a 
dimensionful parameter governing the strength of mixing 
in the Higgs portal. In order to suppress flavor changing 
neutral currents at tree-level, Type-H 2HDMs involve a 
Z 2 symmetry under which Hd —> —Hd and Hu —> Hu- 
We have assumed that this symmetry is softly broken by 
dimensionful couplings, such as by the Higgs portal in¬ 
teraction involving Ba in Eq. 9, and similar terms within 
the 2HDM scalar potential. For simplicity, we assume 
that CP is conserved in the full potential of Eq. 9, which 
implies that Ba is real and that gq and Aq do not develop 


Throughout this work, we will consider values of ma ~ 
100 GeV, niA ~ mn ~ ~ 1 TeV, and 9 ~ 10“^. 

These choices of parameters uniquely determine \Ba\ ~ 
0(0.1) GeV. Therefore, we will be working in the limit 
in which mixing is induced by small off-diagonal terms 
and rriao ^ itt-Ao- ^ result, the light pseudoscalar, a, 
is mostly singlet-like and the much heavier A is mostly 
2HDM-like. 

Stringent constraints on new scalars and pseudoscalars 
can be derived from the results of searches for MSSM 
Higgs bosons at colliders [37-39]. In particular, if a 
has some sizable branching ratio to SM fermions, then 
the production of a a in association with a 6-jet can 
produce distinctive bbb and 6tt events. In our case, 
however, the a has suppressed couplings to quarks and 
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leptons and very large couplings to dark matter, en¬ 
abling collider searches for invisibly decaying light scalars 
and pseudoscalars to provide much stronger bounds [40- 
42]. Even these searches, however, yield extremely weak 
bounds for 9 ^ 10“^ since the production of a de¬ 
pends on its suppressed couplings to SM quarks. In 
fact, for nia ~ 100 GeV in the large tan,0 limit, the 
most stringent constraint comes from the contribution 
to Bs —>■ , which results in the approximate upper 

bound 0 < 0.1 [36]. 

Due to the small mass splitting between xi s-iid X 2 ) 
both of these states can play an important role in de¬ 
termining the thermal relic abundance of dark matter in 
this model. Although we use the publicly available pro¬ 
gram micrOMEGAs [43] to calculate the relic abundance 
numerically, it is illustrative to consider analytic forms 
for the relevant cross sections in the low velocity limit 
(see Fig. 2): 


o-i’(XiXi //) 
o-w(x2X2 //) 
o-u(xiX 2 //) 


ric 

27r 


2nc 

F 

m) 

■ a?iA;^tox ■ 



4to2 — to2 

2nc 

/l^ 

m) 

■ A22A/TOX ■ 



4to2 — to2 


mj^J 


(13) 


3/2 


A5'2Ajmx 

4m2 


where ric = 3(1) for annihilation into quarks (leptons), 
and the widths of a and h should be included when near 
resonance. The first two of these processes are medi¬ 
ated by the exchange of the light pseudoscalar, a, which 
couples to the SM through mixing with the heavier pseu¬ 
doscalar of the 2HDM, XJ = — sin 6 * to/ cot/3/u for up- 
type fermions and XJ = — sin dnif tan f3/v for down-type 
fermions. The last of these processes is mediated by the 
SM-like scalar Higgs boson. For reasons that are similar 
to those described in the previous subsection for the pro¬ 
cess XiXi XiX 2 ) this process is s-wave and contributes 
in the low velocity limit. Here, X 12 is the Xi ~ X 2 — h 
coupling (corresponding to the term A(* 2 ^XiX 2 in the La- 
grangian of Eq. Al), and A^ = —nif/v is the coupling 
between the light scalar Higgs boson and SM fermions. 

The process of thermal freeze-out in this model de¬ 
pends on the hierarchy of the annihilation and coannihi¬ 
lation cross sections described in Eq. 13. If crv(xiXi 
ff), ctu(x 2 X 2 //) > o-n(xiX 2 ^ //), for example, 
each of the two species freeze-out largely independently 
of one another, followed by the decay X 2 —>■ Xi7) which 
increases the final abundance of the xi population. Alter¬ 
natively, if cru(xiX 2 —> ff) is not negligible, these coan¬ 
nihilations will deplete the abundances of both species, 
and the total resulting dark matter abundance. The dark 
matter’s annihilation cross section in the universe today 
can vary significantly depending on which of these pro¬ 
cesses dominates. In the former case, we expect a com¬ 
paratively large cross section, av ~ (4 — 6 ) x 10“^® cm®/s, 
which is in tension with constraints from gamma-ray ob¬ 
servations of dwarf spheroidal galaxies (especially in the 


Field 

Charge 

Spin 

Si 

(1,1,0) 

1/2 

S 2 

(1,1,0) 

1/2 

Di 

(1,2,- 1 / 2 ) 

1/2 

D 2 

(1,2,+ 1 / 2 ) 

1/2 

ao 

(1,1,0) 

0 

Hd 

(1,2,+ 1 / 2 ) 

0 


(1,2,+ 1 / 2 ) 

0 


TABLE I. The field content of the model described in this 
paper. The charges correspond to SU (3)c x SU{2)w x U{1)y. 


case in which 2 to^ is near resonance, but slightly greater 
than TOa, for which the low velocity annihilation rate is 
further enhanced). If the cross section for X 1 X 2 coannihi¬ 
lations is substantial, however, the self-annihilation cross 
section required to generate the appropriate thermal relic 
abundance will be reduced, allowing us to comfortably 
evade this constraint. 

To summarize the major points of this section, if the s- 
channel exchange of the pseudoscalar, a, is to contribute 
to both the up-scattering and the annihilation of the dark 
matter, a must have both large couplings to the dark 
matter and very small couplings to the SM. In our model, 
these latter interactions arise from the mass-mixing of the 
a with the 2HDM pseudoscalar, Aq, allowing the corre¬ 
sponding coupling to be highly suppressed. 


C. Decay and Mass Splitting 

Up to this point, we have assumed that the dark matter 
and its excited state, xi, 2 ) are gauge singlets. While 
this is sufficient to obtain the desired rates for both up- 
scattering and annihilation, we must also require that 
the excited state, X 2 , decays with a lifetime that is much 
shorter than the age of the universe. This requires the 
dark sector to couple to a charged state appearing in the 
loop-diagram responsible for the decay X 2 —> Xi7- 

A cosmologically short lifetime for X 2 can be easily 
accommodated by mixing a small SU(2)w doublet com¬ 
ponent into the dark sector, allowing xi ,2 to couple to 
the Higgs doublets, Hd,u, (and their associated H^) as 
well as to . Additionally, an active degree of freedom 
that is allowed to mix with xi ,2 generically introduces 
new charged states into the dark matter sector, analo¬ 
gous to charginos in the MSSM, that can enter into the 
loop-induced decay as well. The introduction of a small 
doublet component into the dark sector can also provide 
a natural explanation for the 3.55 keV mass splitting be¬ 
tween the states xi and X 2 - 

In this regard, we follow closely the approach laid out 
in Ref. [44], wherein a vector-like pair of 2-component 
Weyl fermion SM gauge singlets, and S 2 , and a vector¬ 
like pair of Weyl fermion SU{2)w doublets, Di and D 2 , 
are introduced, the latter of which are assigned hyper¬ 
charge t 1/2. Unlike in Ref. [44], however, which only 
considers long-lived primordial decays with mixing in- 
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FIG. 3. Feynman diagrams for the decay X 2 —>■ Xi7- Simi¬ 
lar diagrams in which the photon is emitted off the x^ S'lso 
contribute. 


troduced by interactions with the SM Higgs donblet, we 
will consider interactions involving the two Higgs dou¬ 
blets and much shorter lifetimes for the excited state. 

In addition to the bare mass terms, in general, the 
singlet and doublet degrees of freedom can couple via 
Yukawa terms to one or both of the Higgs doublets, 
Hd,u- We will assume that the dark matter sector re¬ 
spects the Z2 symmetry of the 2 HDM potential, which 
we enlarge to include S'i_2 —t —S'!,2 (in addition to the 
usual Hd u —t THd,u)- As a result, the dark matter can 
only directly couple to Hd- We summarize our model’s 
particle content in Table I. The dark sector Lagrangian 
contains the following terms: 

—C D MsSiS2 + MJJD1D2 + iy30,08182 + iyDOoDiD2 
+ yii8iDiHd + y 2 i 82 DiHd 
+ y22 82H^d ■ ^2 + yi28iH^ ■ D2 + h.c., ( 14 ) 

where 2 -component Weyl and 8U{2)w indices are im¬ 
plied. For simplicity, we take all of the couplings in Eq. 14 
to be real and introduce an additional Z2 symmetry on 
the dark sector such that the lightest fermionic state is 
a stable dark matter candidate. For model building in a 
similar direction, see Ref. [ 45 ]. 

The dark matter 8U{2)w doublets are parametrized 
as: 



where V12 and E12 are the neutral and charged compo¬ 
nents of the doublets, respectively. Ei and E2 will mix to 
form an electrically charged Dirac fermion of mass Mo, 
which we will label as x^- Since chargino searches at 
LEP generally exlude new charged fermions lighter than 
100 GeV, we will take > 100 GeV. After elecroweak 


symmetry breaking, the neutral degrees of freedom, S'i_2 
and 1^1,2) will mix according to the following mass matrix 
(in the 81-82-1^1-1^2 basis): 


/ 0 Ms ^yiiVd ^yi2Vd\ 

Ms 0 ^y2iVd ^y22Vd 

-^yiivd ■^y2ivd 0 Mo 

\^yi2Vd -^y22Vd Mo 0 y 


( 16 ) 


The lightest mass eigenstate of Mq will be the stable 
dark matter candidate, xi, and the second lightest mass 
eigenstate the excited state, X2- The gauge composition 
of Xi,2 can be written as: 


Xi,2 — ^sfSi -\- Ng^82 + Nlfvi NI^V2 ■ ( 17 ) 

A large doublet component of xi,2 is severely restricted 
by direct detection experiments because this usually in¬ 
troduces large couplings to the SM-like Higgs and to the 
Z. We therefore focus on the case in which xi and X2 are 
largely singlet-like and « Ms ^ Mo- In this 

small mixing limit, it is possible to derive approximate 
forms for the mixing angles: 


N. 


1.2 

Si 


T 






1 

'71 


( 18 ) 


Vd 


■ [ (2/11 ± 2/12) Mo ± (2/21 ± 2/22) Ms ] 


~ 2 (M^^- Af|) I- ^ ^ ^ 

Eurthermore, in this limit, the mass splitting, Srriy, = 
^X2 ~ he approximated as: 


6m^ 


Ml-Ml 


(2/112/12 + 2/212/22) A 7 £) 

+ (2/112/21 +2/122/22) Als 


( 19 ) 


Erom Eq. 19 one can see that if either yu = 1/22 = 0 
or yi2 = 2/21 = 0 , then Sm^ = 0 . This can be under¬ 
stood from the symmetries of the Lagrangian of Eq. 14 
as follows. The kinetic terms of Ai, 82, Di, D2 possess a 
27 ( 1 )"^ symmetry. This is broken down to U{l)s x U{l)o 
by the bare masses Ms and Mo of Eq. 14 . At this point, 
27(1)5 X 27(1)51 guarantees that the mass eigenstates of 
Mq in Eq. 16 will split into two separate degenerate pairs. 
Once the Yukawa couplings yu, yi2, 1/21, 2/22 are turned 
on, both of these 27 ( 1 )’s are further broken. However, in 
the limit that either yu = 1/22 = 0 or 1/12 = 1/21 = 0 holds, 
then 27(1)5 X 27(1)5) is restored, and the mass spectrum 
once again decouples to two pairs of mass degenerate 
states. This argument holds to all orders in perturbation 
theory. 

To obtain a mass splitting as small as 3.55 keV, 
Yukawa couplings on the order of ~ 10 “^ x 

(tan/ 3 / 50 )(M 5 )/ 700 GeV)^/^ are generally required. As 
we will see in the next section, however, order one values 
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FIG. 4. The lifetime for the decay X 2 —^ Xi7: for Sm^ = 3.55 
keV, ^ MD,'mu±, mjj± = 1 TeV, tan/? ^ 1, yn « 
J/21 ~ 2/22 ~ —J/12 ~ 2.5. 


for these quantities are necessary if X 1 X 2 coannihilations 
are to be efficient enough to obtain the desired thermal 
relic abundance without conflicting with constraints from 
gamma-ray observations of dwarf galaxies. This tension 
can be resolved if there is a cancellation in Eq. 19, such 
as arises for the choice yu ~ 1/22 — 2/21 — —yi 2 , with a 
small level of non-degeneracy needed for an ©(keV) mass 
splitting. This is technically natural and can arise from 
symmetries in the ultraviolet, which are expected to be 
broken at the two-loop level (for a similar approach, see 
Ref. [44]). 

After electroweak symmetry breaking, xij X 2 , and x^ 
couple to the charged bosons and W^, and a mass 
splitting, 6m^, is induced. As a result X 2 can decay via 
a radiative 2-body process X 2 Xi7 where x^ and the 
charged bosons are exchanged in the loop (see Fig. 3). 
In the simple limit that yu = j /22 = J/ 2 i = —yi 2 = y, 
mH±,MD m^, and tan/3 ;g> 1, the loop con¬ 
tributes negligibly, and this decay width is given by: 


r 


e'^y^ 

25671® (m^± — 


1 - 


mjj± - Ml 


In 


Ml ) 


2 


( 20 ) 

A more general expression for this width is given in Ap¬ 
pendix B. We find that the most viable parameter space 
leads to lifetimes for X 2 that are on the order of an hour 
(see Fig. 4). Furthermore, the 3-body decay to neutrinos, 
X 2 remains kinematically open via an off-shell 

Z boson. This later process, however, remains subdom¬ 
inant since it is suppressed by two additional powers of 
5m^ [44]. 


In the region of parameter space under consideration, 
our model is not significantly constrained by observa¬ 
tions of the cosmic microwave background (CMB) or of 
the light element abundances. In particular, for X 2 life¬ 
times less than 10^^ s, the most stringent bounds from 
the CMB are on the XiXi annihilation cross section, and 
are weaker than the values considered in our study [46]. 



FIG. 5. An example of the parameter space of our model, 
with |AJi, 22 | = 6, i/d = 0, yu = 2/21 = 2/22 = 2.5, 7/12 « -2.5, 
5m^ = 3.55 keV, 0 = 3 x 10“®, tan/? = 50, nia = 120 GeV, 
and mA ~ rnu = mH± ~ 1 TeV. The region yielding a 
thermal relic abundance equal to the cosmological dark mat¬ 
ter density = 0.1199 ± 0.0027 (solid black line) passes 
through the region that can generate the observed 3.55 keV 
signal (labelled “XDM”) for values of Ms « 61 GeV (about 1 
GeV above the a resonance). Also shown is the region capable 
of generating the Galactic Center gamma-ray excess (labelled 
“GCE”, corresponding to au = 5 x 10“^^ to 5 x 10“^® cm®/s) 
and the regions that are excluded by gamma-ray observations 
of dwarf galaxies [49] or by LUX [50]. 


Constraints from Big Bang Nucleosynthesis (BBN) are 
generally expressed in terms of the rescaled electromag¬ 
netic energy released, as a function of the X 2 lifetime. 
For lifetimes less than 10® s, the upper bound is derived 
from the measured deuterium abundance [47, 48]. Our 
model, however, predicts values for this quantity that 
are many orders of magnitude smaller than the existing 
upper limit. 


IV. RESULTS 


Due to the sizable number of free parameters in this 
model {Ms, Md, yii, 2 / 22 , 2 / 12 , 2 / 21 , yo, ys, a, P, Wao, 
iriAo, Ba), a wide range of phenomenology can emerge. 
We will simplify this to some extent by focusing on the 
parameters which yield ~ 3.55 keV and ~ 60 
GeV, which are within the range capable of generating 
both the 3.55 keV line and the Galactic Center gamma- 
ray excess [29]. 

In Fig. 5, we show an example of a slice of the pa- 


















rameter space in this model. Here, we have adopted 
1 ^ 11 . 22 ! = ^, VD = 0, yii = y2i = 2/22 = 2.5, yi 2 ~ -2.5 
(with 2/12 fixed to give Srriy^ = 3.55 keV throughout the 
Md — Ms plane shown), 9 = 3 x 10“^, tan/3 = 50, 
Too = 120 GeV, and uia = rnn = rnH± = 1 TeV. The 
region yielding a thermal relic abundance equal to the 
cosmological dark matter density (solid black line) passes 
through the region that can generate the observed 3.55 
keV line (labelled “XDM”) for values of Ms « 61 GeV 
(about 1 GeV above the a resonance). In particular, to 
provide an adequate fit to the 3.55 keV signal we demand 
that the up-scattering rate fits the measured fluxes within 
the < 3 contour of Ref. [17] for promptly decaying 
excited states. Also shown is the region that is capa¬ 
ble of generating the Galactic Center gamma-ray excess 
{av = 5 X 10“^^ to 5 X 10“^® cm®/s), as well as the 
regions that are excluded by gamma-ray observations of 
dwarf galaxies [49] or by direct detection experiments [50] 
(see Appendix C). Constraints from the invisible width 
of the Higgs do not restrict any of the parameter space 
shown. 

The process of X 1 X 2 coannihilation plays an important 
role in the determination of the thermal relic abundance 
in this region of parameter space. In particular, in the 
limit at hand, ^ 12/^11 ~ ^Md/Ms ^ 1, large coannihi¬ 
lation rates are possible without large elastic scattering 
cross sections with nuclei (for details, see Appendices A 
and C). In addition, resonant annihilation and coannihi¬ 
lation through the a and the SM-like Higgs significantly 
deplete the thermal abundance for Ms — 60 GeV and 
62.5 GeV, respectively. The relic abundance is relatively 
insensitive to the parameters scanned over in Fig. 5, and 
is within approximately an order of magnitude of the 
measured quantity over the entire plane shown. Lastly, 
we note that up-scattering and decay rates are insensitive 
to the parameter 9. If we had chosen to set this quan¬ 
tity to zero (decoupling a from the pseudoscalar of the 
2HDM), we can still generate the 3.55 keV line, but with¬ 
out a mechanism to produce the Galactic Genter gamma- 
ray excess and without constraints from gamma-ray ob¬ 
servations of dwarf galaxies. 


V. SUMMARY AND CONCLUSIONS 

It has been previously proposed that dark matter scat¬ 
tering into an excited state (eXciting Dark Matter, or 
XDM) could be responsible for the 3.55 keV line ob¬ 
served from Galaxy Clusters without conflicting with the 
lack of such a signal from dwarf galaxies [15]. Such a 
model could also potentially generate Fermi’s gamma- 
ray excess from the Galactic Center. Most of the XDM 
model building discussed in the literature has focused 
on scenarios in which the dark matter interacts through 
a light mediator, with no significant couplings between 
the dark sector and the Standard Model. Here, instead 
of hidden sector, we have considered a model in which 
the dark matter directly annihilates into Standard Model 


fermions through the near resonant exchange of a pseu¬ 
doscalar, a, which also efficiently mediates the process of 
up-scattering, XiXi —>■ X 2 X 2 - This pseudoscalar is a mix¬ 
ture of a Standard Model singlet and the pseudoscalar 
appearing from a two-Higgs doublet model. The dark 
matter itself is a mixture of two Standard Model gauge 
singlets and the neutral components of two SU (2)iy dou¬ 
blets. This allows us to generate a 3.55 keV mass splitting 
between the two lightest mass eigenstates, and enables 
for the rapid decay of X 2 —>■ Xi7- 


We have identified regions of parameter space in our 
model that can simultaneously generate the 3.55 keV line 
and the Galactic Genter gamma-ray excess, while remain¬ 
ing consistent with all constraints from colliders, direct 
detection experiments, and gamma-ray observations of 
dwarf galaxies. Coannihilations between xi and X 2 can 
play an important role in determining the thermal relic 
abundance of dark matter in this model. 
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Appendix A: Higgs and Gauge Couplings 


In this appendix, we provide analytic forms for the couplings of XiXi> X 2 X 27 and X 1 X 2 to the light Higgs bosons (a, 
h) and to the Z. These couplings are defined according to the following terms in the Lagrangian: 

C D Kiaxiil^Xi + A22aX2*7^X2 + K2^Xiil^X2 + >^\ihxiXi + A22^X2X2 + A^2/iXiX2 

+ 9llZ^,Xll^l^Xl + ff22^MX27^7^X2 + 912 Z^Xil^ 1 ^X 2 - (Al) 

As discussed in the text, the field X 2 requires a transformation of the form X 2 * 7 ^X 2 in order to ensure a positive 
value for its mass term. After this field redefinition, the above Lagrangian appears as follows: 

£ D Xhaxiil^Xi - A22aX2*7^X2 - A!j*2aXiX2 + Anft-xiXi - A22^X2X2 + Ai2/iXi*7®X2 
+ ffii-^AiXi7^7^Xi + 922 Zf,X 2 l^ 1 ^X 2 + i9i2Z^Xil^X2- (A2) 

The couplings are given by: 


= cose[ysNlNl + yoNl^Nl^] + + y2iNl^Nl^ - y22NlNl^ - 2/i2iV^,<] 

A^2 = COS0 [j/sA^|,A^i + VdNIX^] + [yiiNlNl^ + y^iNlNl^ - y22NlNl^ - yi2NlNl^] 

A?2 = cosB[ys{NlNl + + y^(A,\A^^ + A^^A^)] + 

+ y2i{NlNl^ + Aj^ A^J - y22{NlNl^ + AJ^ A^J - 2/12 (A^^A^^ + A|^ A^J] 
h _ cos/3 

h _ cos/3 

h _ cos/3 


+ y2iNl^Nl^ + y22Nl^Nl^ + yi2Nl^Nl^ 
+ y2iNlN^^ + y22NlNl^ + yi2NlNl 


911 = 

922 = 

912 = 


9 

Acw 

9 

Acw 

9 

2 cw 


2 /ii(A^Si< + NlKj + // 2 i(A 1 ^A, 
(Aif-(A4f 


+ y22{NlNl + A|^ Ai) + yi 2 (iV^,< + iV |,0 


(a 2 j 2 _ (iv 2 )2 


A^ - A^ A^ 

Vi V 2 ^ V 2 


where the mixing angles are defined in Eq. 17, and g and cw are the SU{2)w coupling constant and cosine of the 
Weinberg angle, respectively. Note that in the limit of small Oq-Ao mixing, A“j 22 ~ ±j/sCOS0/2. 


Appendix B: Decay 

In this appendix, we provide formulae describing the decay X 2 Xi7- We note that a convenient gauge choice is 
non-linear R gauge due to the vanishing of the ylE+G” vertex (see e.g. Ref. [51]). The width for this process is given 
by: 


r(x 2 


X 17 ) = 


Stt V Tn 


A 


''X 2 


(Bl) 


where A = . The effective coupling in this expression is given by: 


^efT — 


eei 

Stt^ 


- (eie2ffiL52L - 9iR92R)[c2m^^{Ig2 - Jg- Kg) + tim^^{Jg - Kg)] - 2MD{eigiLg2R - C2giRg2L)Jg 


+ ^(AilA2l — eie2AiflA2fl)[e2TO;^,2(/s2 — Kg) — eim^^Ks] + -M£)(e2AiLA2K — eiAi_RA2L)/s 
+ ^(AfiA^^ — eie2AfflA^fl)[e2TO;^2(/g2 — Kg) — eimy^^Kg] + - Afp(e2Af^A^^ — eiAf^A^^)/g 


(B2) 
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where ei and 62 are the signs of the first two eigenvalues of the mass matrix, Mq, as given in Eq. 16. The couplings 
in the above expression are given by: 

= giR = g 2 L = --^Nl^, g 2 R = -^Nl^, (B3) 

Ail = -sin/3(?/22^S2 +J/i 2 ^Si)> ^ 2 L = - sin/? (j/ 22 iV|^ +yi 2 ?V|J, 

\iR = -sin/3(2/ii7V^^ +J/ 2 i^lJ, A 2 i^ = - sin/3 +921^%^), 

>'IL = cos /3 (2/22^^52 + J/l2^Si)> ^2L = COS (3 {y22N^^ + yi2N^J, 


Afi? = cos/3 ( 2 / 11 + 2 / 2 iiV^J, = cos/3 ( 2 / 11 3V|^ + 2/2i3V|J. 


Lastly, Eq. B2 contains a number of integrals, defined 
as follows: 


A 


42 = ^ 


4 - A 


■logATg, 


da; log Xg , 


logX, 


atomic mass A is given by: 


Zfp + (A1 - Z)/, 


where is the reduced mass of the system and the 

nucleon level couplings are given by: 

f _ ™ r(p-”) I A f (p>") _4_ 

Jp,n-mp,n 2^ St, + 27 Atg 2^ > 

^ q—u,d^s ^ q—c,b,t ^ 


A 


1,2 = -r 


A 


■ logA:^, 


da; log Xs 


■logAf', 


dslg = - ^(1 + M|)/g + m^Jg - m'^^Ig2), 


Kg = ~ ^ (1 Af^Ig + ?Ti^± Js — 


mlJs2) 


where 


_ Ml,x + m'^{l - x) - - x) 

® M|,a; + m^(l — x) — TO2^a;(l — x) ’ 

^ + Ad|,(l - x) - to2^x( 1 - x) 

® m^x + M|)(l — x) — to2^x( 1 — x) ’ 
^ _ MjgX + rn\j± (1 - x) - rri^^^xil - x) 
* M|,x + (1 — x) — m^^x(l — x) 

^ m^±a; + M|)(l - x) - to2^x( 1 - x) 
TO^±x + M|,(l — x) — m2^x(l — x) 


1 All I Af^g/3 

2 2 ’ 

TUq V [ 

where qp = cot/? (—tan/?) for up-type (down-type) 
quarks. 

In addition, direct detection experiments can also de¬ 
tect inelastic events, xiA —>• X 2 ?V- For Snix ^ x^/a;^,Ar/2, 
the dark matter particles typically possess enough kinetic 
energy to scatter into the excited state. The cross section 
for inelastic scattering is given by: 


F Zfp + {A-Z)U 


where in this case, 


gsin 6 wgi 2 


(cot div — 3 tan 9^) 


gs\n9wgi2 , ,, . ^ 

Jn = --(cot 9w + tan Ow), 

and the following factor accounts for the kinematic sup¬ 
pression associated with inelastic scattering: 


- 2(to2^ -I- m%)s + 

— 2(m^, + m4)s -I- (m^^ — _ 


, (C6) 


Appendix C: Direct Detection 

The elastic scattering cross section between the dark 
matter, xij and a nucleus with atomic number Z and 


\/s « niN 1 -I- 
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where v ~ 10“^. Note that _F —)• 0 in the limit of f —> 
riinin = below which inelastic scattering 

is not possible. In contrast, for the mass splitting and 


masses under consideration in this paper (and for typical 
dark matter velocities in the local Milky Way), v ^ rimin 
and F ~ 1. 
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